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Analytic phenomenological scaling is carried out for the random field Ising model in general 
dimensions d using a bar geometry. Domain wall configurations and their decorated profiles and 
associated wandering and other exponents (C,7, 5, fi) are obtained by free energy minimization. 
Scaling between different bar widths provides the renormalization group (RG) transformation. Its 
consequences are (i) criticality at/i = T = 0ind<2 with correlation length £(h,T) diverging 
like £(h, 0) oc h- 2/(2 - d) for d < 2 and 0) oc exp[l/(ci7/i 7 )] for d — 2, where ci is a decoration 
constant; (u) criticality in d = 2+e dimensions at T = 0, h* = (e/2ci) 1 ^" f , where £ oc [(a— s*)/s]~ 2t ^ 1 , 
s = h 1 . Finite temperature generalizations are outlined. Numerical transfer matrix calculations 
and results from a ground state algorithm adapted for strips in d — 2 confirm the ingredients which 
provide the RG description. 



*Electronic address: dmoore@thphys.ox.ac.uk 
"^Electronic address: sldq@portela.if.uff.br 



1 



The random field Ising model (RFIM) , which is closely 
related to dilute antiferromagnetic Ising systems in a uni- 
form field [Q, has provided a sequence of challenges to 
theoreticians as well as experimentalists [||-|| . The even- 
tual resolution of the puzzle concerning its lower crit- 
ical dimension confirmed the validity of domain wall pic- 
tures 0. Attempts to elucidate the critical behaviour 
at the lower critical dimension, d = 2, by normally very 
powerful numerical techniques J^— [To|] have encountered 
difficulties related to the apparently anomalously severe 
divergence of the correlation length there. The phase di- 
agram boundaries are not known for general dimensions, 
and the critical behaviour in three dimensions HJl3| is 
presently not understood. 

Yet methods for describing individual domain walls 
have been provided fl3|- 16f| , and their exploitation in 
standard geometries fll5|l6fl has been very effective, 
within the limits of a scheme without renormalization. 

The aim of the present work is to use domain wall 
pictures in a bar geometry to build up, through finite 
size scaling [ [L7|]l8| , a renormalization group (RG) de- 
scription from which the critical properties can be ob- 
tained. For strip geometries in d = 2, we also present 
preliminary free energy results from a numerical trans- 
fer matrix procedure Jl9| from which quantities such 
as the free energy and correlation length can be esti- 
mated. Though presently short of providing a full self- 
contained numerical phenomenological RG description, 
these together with and a direct numerical solution for 



the ground state |20| allow a check of the analytical 
domain wall arguments from which the renormalization 
scheme is built up. 

The method exploits RG transformations obtained by 
finite size scaling in the bar geometry shown in Fig. [I]. 

We begin with the case of zero temperature T, which 
most easily illustrates the procedure and already includes 
several important results and considerations. Here the 
basic quantity is the energy U per site. Throughout the 
letter the no-wall term is subtracted, and hypercubic or 
continuum cases are considered. For flat domain walls 
with average spacing £ L , U becomes U F = WL 1 ^ d /£ L 
where 



W = J{2L d ~ 1 - Acoh^LL^ 1 )^ 2 ) 



(1) 



This is made up of a surface term and one proportional 
to (volume) 1 / 2 . The variable h is the rms random field 
in units of the exchange interaction J, d is the dimen- 
sionality and Co is a constant of order 1. Minimizing U F 
with respect to £l gives the equilibrium values £ F , U F : 



{c h)~ 



U F = -2(c h) 2 /L 



d-l 



(2) 



As discussed later, these primitive results are asymptot- 
ically correct for small h, where (i, > L, and are con- 
firmed by comparison with the numerical transfer matrix 
(TM) results in the appropriate limit. 



The RG transformation of parameters (here h — > h' = 
Rb(h)) under rescaling of L by b is obtained from the 
phenomenological scaling relation [l7|JT^ | 



£ L (h')/L = {[ bL (h)/(bL) 



Using £ F , as in Eq. |[ gives 



ti = R F (h) EE b^^h 



(3) 



(4) 



This result has been previously obtained in the field- 
theoretic scheme of Grinstein and Ma pjj ; see also Ref. 
p2| . Eq. |i| has an unstable fixed point at h = for 
d < 2. Here the eigenvalue 6( 2_d )/ 2 implies that the bulk 
correlation length diverges for small /iatT = 0ind<2 
like 



£ oc h 



-2/(2-d) 



(5) 



The transformation (|J) is marginal for d — 2. So, 
subdominant terms are required, even at small h, T = 0, 
for d — 2, to deal with the marginality. These come 
from "decoration" perturbations of the flat walls [O [Tg| . 
These can be built up from smooth shape modifications 
of base length a and height 6(a), superimposed for all 
scales a between L and a cutoff value determined by the 
lattice constant (=1), see Fig. |^ . The minimum energy 
s(a) and minimizing height b for a smooth decoration of 
base scale a are of the form 



e(a) oc -h 7 a s , 6(a) oc Wa^ 



(G) 



where C is a wandering exponent (lj^l|,g3|,gj| , and the 
exponent 7 will be particularly significant below. The 
minimization gives (general d): 



4 d+1 2 5- 

M,m,C) = (3.— ,3, -3- 



or (2, 1,2, 3-d) , 
(7) 



where the first result applies for the continuum case if 
b < a, and the second one applies to the continuum if 
6 > a, and to the lattice case. So, in the continuum case, 
if a is large the first result applies for d > 2 (largest 5, 
i.e. lowest energy) and the second for d < 2. 

Since in Eq. ||e(a) is negative, decorations on decora- 
tions occur down to a smallest base length a m such that 
min(a m ,6(a m )) « 1. The resulting contribution of all 
decorations to the energy of a single wall of base length 
L is, for d — 2, 



-Jci/i 7 Lln(L/^c 2 ) 



(8) 



where C\ and c 2 are constants of order 1. The effect of 
decorations is therefore to replace the wall energy W by 
W + X. Minimizing the energy per site with respect to 
wall separation then gives equilibrium values U D , £^ for 
the decorated wall generalization, and hence (via £f / ' L) 
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the generalized phenomenological scaling equation. For 
d near 2 (the marginal case), this is: 

ti = W/ 2 - d)/2 [l + ci/i 7 ln6] . (9) 

From Eq. 0, 7 = 4/3, 2 for continuum and lattice cases 
respectively. The transformation given by Eq. ^ has the 
proper semigroup character, in its range of validity {h 
small, d near 2, T = 0). As will be seen, this range 
includes the critical effects for d at or just above 2. 

For d = 2, Eq. || shows that the h = fixed point 
is marginally unstable, as expected at the lower critical 
dimension. Analyzing the equation by standard RG pro- 
cedures gives the two-dimensional bulk correlation length 
diverging for small h and T = as 

In 2 + e dimensions, e small and positive, the RG transfor- 
mation Eq. |^ yields an unstable zero temperature fixed 
point at small non-zero field h* , as well as the associated 
divergence of £ : 

h* = {e/2 Cl ) l h , £ cx [(s - S *)/,sr 2 ^, 8 = K< . (11) 

This divergence crosses over to the behaviour (|l^) as 
e — > 0. In Eqs. ([l0]) and ([ll]) the exponent 7 differs be- 
tween lattice and continuum systems. This feature has 
previously been noted in the interface context |lj] . 

For larger fields, the flat wall description requires more 
drastic generalizations. Increasing h at fixed L would, 
according to Eq. 0, eventually give fx < L. So, by 
then the effects of domain boundaries in all directions 
should have been allowed for. This generalization ]25| ] 
yields a complicated form which confirms Eq. |l| and its 
consequences in the low field regime where fx 3> L, but 
which for (i«Jj gives 

U = J[c 1 &-c 2 h$' t '- 1) ] , (12) 

where c±, C2 are geometric constants and <fi crosses over 
between 1/2 and 1 as fx approaches order 1. For d = 2, 
T = the result of allowing for this, and the decoration 
effects, is as follows. As h increases the low field result for 
U arising from Eq. || (namely — 2(coh) 2 / L) picks up a fac- 
tor A(h, L) = [1 + c\ K 1 ln(L/i M C2)], then goes through a 
complicated intermediate regime, and finally for h larger 
than about 1 becomes linear in h and independent of L. 

A check of the basic (flat wall) ingredients in the do- 
main scaling description is provided by comparing these 
results for U with those provided by a numerical trans- 
fer matrix calculation for the d = 2 RFIM |lS| ], Here 
the configurationally averaged free energy is enumerated 
exactly for very long (~ 10 5 lattice parameters) strips 
of finite width L. Fig. ^ gives the field dependence of 
the negative of the free energy per spin /, for relatively 
small widths L = 2, . . . 7, at temperature T/J = 0.1 



which is sufficiently low that / is essentially U . The 
main characteristics of the domain wall results for —U 
{h 2 /L crossing over to h/L°), may be seen to be present. 
Indeed the analytical results, taking cq = 0.4614 in Eq. 
|2|, fit the numerical data for L = 4 — 9 within the very 
small error estimated from the data fluctuations in the 
low field regime h < 0.5 for T/J = 0.1, and also for higher 
temperatures using extensions to the domain wall theory 
outlined immediately below. However, the fit is not suf- 
ficiently sensitive to the decoration terms to confirm the 
theoretical value 7 = 2. Other comparisons discussed 
later show more direct effects of the decoration terms. 

The finite temperature generalization involves entropic 
terms. At the most primitive level (flat walls only) the 
entropy per site is: 

S F = -k B L^ d \x\nx + (1 -a?)ln(l - x)], x = 

(13) 

coming from the number of ways of selecting Cx out of 
C sites for the placement of the walls. Adding —TS F to 
the flat wall internal energy U F and minimizing gives the 
equilibrium flat wall separation £ F and the free energy 
f F per site. £ F satisfies 

c h(Z F /L d - 1 ) 1 / 2 + (L 1 - d /2K)]n(Z F -l) = l , (14) 

where K = J/k B T. 

This result is valid if fx, 3> L, which is the case if 
both hlP-®/* and (l/K)V- d hxL are small, sufficient 
to encompass Eqs. |^, ^ and [5] and their low-temperature 
generalizations. Otherwise domain walls in all directions 
are needed, and they increase the total entropy, asymp- 
totically to C3^ d ln(c4fi) where 03,4 are constants. And 
in d = 2 decoration terms are required to resolve the 
marginality at low T. They affect both the energy as dis- 
cussed above (taking K in Eq. [l4| to KA(h, L)), and the 
entropy. Their entropy contribution for each vertical wall 
of base length L in d = 2 is of the form c^Lh^ ln(ceL). 
The resulting temperature-dependent generalizations of 
Eqs. |, |, |, [l(] and |TT] will be given elsewhere |j§, 
together with phase boundaries and comparisons be- 
tween analytic and numerical transfer matrix and Monte 
Carlo p6| results for finite temperature free energies. 

However, we briefly present here evidence for the 
marginality-breaking domain wall roughening given in 
equation (Q) and used to derive key results such 
as (§)) (H)j ( |l0| ) and (|ll"l). Specifically, we demonstrate 
that, for a gaussian distribution of random fields, the ex- 
ponent [i approaches 2 as L grows large. The max-flow 
algorithm of Ogielski [^o| was implemented on a strip 
geometry of dimensions L x 1000 and the ground state 
was generated for 100 independent random field configu- 
rations of fixed standard deviation h. For each spin con- 
figuration, the rms width of the domain walls was mea- 
sured, and this statistic was averaged over all the domain 
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walls of the 100 field configurations. We expect b rms (L) 
to scale like b(a) in equation (^). Figure || justifies this 
prediction for large L. 

Detailed results, discussions of the operational proce- 
dures will be given elsewhere |ll||25|,^6|] . The results are 
here used merely to support the domain scaling approach. 

In conclusion, an analytic phenomenological scaling 
approach has been constructed using domain considera- 
tions. The resulting RG transformations have been used 
to find critical properties, including the correlation length 
behaviour at and near the marginal (lower critical) di- 
mension d — 2. Numerical transfer matrix free energy 
comparisons have confirmed, in d — 2, the basic domain 
picture used in the scaling. And numerical ground state 
wall roughening investigations have quantitatively tested 
a particular theoretical prediction (the exponent /i) for 
the decoration ingredient, which is so crucial for the crit- 
icality in 2 and 2 + e dimensions (see Eqs. |l0|, |TT| ). 

Generalizations of the domain scaling approach to re- 
lated vector model systems and to dynamics are being 
pursued starting from existing domain wall pictures (re- 
spectively Refs. [0 and |||). 
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FIG. 1. Domain walls of separation £ i; in a d-dimensional 
bar of transverse scale L. 

FIG. 2. Domain wall decorations (d — 2 for simplicity) of 
base length a, height b(a): (i) for continuum; (ii) for lattice; 
(Hi) superimposed at successively smaller scale. 

FIG. 3. Free energy for the d = 2 RFIM, from nu- 
merical transfer matrix calculations for long strips of widths 
L = 2, 3, ... 7, as a function of h/J at T/J = 0.1. The be- 
haviour crosses over from h 2 /L to h/L°, as predicted by do- 
main wall arguments and used in scaling. 

FIG. 4. The domain wall roughening exponent fi. This 
has been determined by least-squares fits to semi-log plots of 

brms(L) VS h. 
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